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The inverse scattering method for cyhndrical 
gravitational waves 
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The initial- value problem for cylindrical gravitational waves is stud- 
I led through the development of the inverse scattering method scheme. 

The inverse scattering transform in this case can be viewed as a trans- 
formation of the Cauchy data to the data on the symmetry axis. The 
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Abstract 



qh| Riemann-Hilbert problem, which serves as inverse transformation, is 

^ ' formulated in two different ways. We consider Einstein-Rosen waves 

bX)' to illustrate the method. 

> : 

^ ' 1 Introduction 

In this paper we study the cylindrical gravitational waves model with two 
polarization modes. The space-time of the model possesses two commuting 
one-parameter isometry groups of which one is isomorphic to R and the other 
to SO [2), and is non-stationary solution of the vacuum Einstein equations. 
The generators of these groups are called Killing vector fields. In the case of 
one polarization (both Killing fields are hypersurface orthogonal) this model 
is also known as Einstein-Rosen waves. The theory of Einstein-Rosen waves 
is linear, however, still physically interesting and even in recent time provides 
us with new important results both in classical and quantum scope Ref. 
where one finds additional discussion and reference list. The polarized model 
is non-linear and hence far more difficult to study. 
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In two Killing reduction of vacuum or electro-vacuum general relativity 
Einstein equations belong to the class of integrable equations. By the other 
words they can be written as a compatibility condition of an auxiliary linear 
system. Although in this paper we use the linear system proposed in Ref. 
it is worth mentioning that there exits essentially different approach Ref. [^] . 

The best tool to study integrable equations is the Inverse Scattering 
Method (ISM). There are many works where the ISM is applied for con- 
structing exact solutions and there are a few papers where the physically 
relevant models are considered. The author knows only Ref. where the 
model of colliding gravitational plane waves was investigated and Rev. 
where the solutions with disconnected horizon were analyzed. 

The main goal of the present paper is to develop a framework to study 
the cylindrical gravitational waves model. 



2 Field equations and boundery conditions 

The metric of the space-time with two commuting space-like Killing fields 
can always be chosen locally in the following form 

ds'^ = dsji + Xd(f)'^ + 2Wd(f)dz + Vdz\ (1) 

where 

dsjj = JabdXadXb 

is the metric of an orthogonal to both Killing vectors two-surface. The vari- 
ables X, W, V depend only on the coordinates Xa adapted to this surface. 
Let 

f V W \ 2 , , 

9=[w X J' ^ ='^^^^- 

Then, the first group of Einstein equations can be written as 

d*pdgg^^ = 0. (2) 

Here * is the Hodge operator with respect to the metric dsjj. The function p 
has a space-like gradient hence it can be used as a radial coordinate. Let the 
time coordinate be dual to space one, dt = *dp. Coordinate chart constructed 
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is often called the Weyl canonical coordinates. We employ it throughout the 
paper. In these coordinates the 2-metric ds'jj is given by 

dsh = f{p, t) (-de + dp') , f{p, t) = (3) 

The last formula is a definition of the function 7. Then the system takes 
the form 

- {pg,tg'^),t + {pg,pg'^),p = 0. (4) 

We impose the following boundary conditions on g 



{ V p'W ^ 



where V, W and X are smooth functions not equal to zero and 

V = l + J/p + 0(l/p2), X = 1 - J/p + 0(l/p2), W = 0(l/p3) 

at p = 00. Here J is a constant independent on t. 

An alternative formulation of Einstein equations (Q) exists. For this, let 
us introduce matrix potential 

dH = p*dgg-\ Y = Hu. (6) 

Here Y is the Ernst potential that is determined by the transition Killing 
vector field. The functions Y and V are independent dynamical variables 
which set a solution of (0) completely. Potentials Y and V satisfy the Ernst 
equation which we will not use in this paper. 
Using d^) one can easily prove that 

-1 ( dtY\ _i 



at p = and 



oa a-' f \ _ / 0(l/p) 0(l/p3) \ 
^^'^^ " 1^ 2 j - 1^ 0(l/p) 0(l/p) ) ' 

,-1 _ 0(l/p) 0(l/p3)\ 

- I o(i/p) o(i/p) j ' 
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at p = oo. 

The second group of Einstein equations allows one to determine the co- 
efficient 7 from the matrix g. From (^) it is possible to show that = as 
p = 0. The solution is regular at the symmetry axis iff 7 = for p = 0. How- 
ever, in this case 7— s>7o7^0asp— »oo and the conical singularity appears 
at the space infinity. Its angle of deficit can be treated as the energy of the 
system [||. It is interesting to note that the deficit angle of the conical sin- 
gularity at the symmetry axis of the stationary solutions with disconnected 
horizon has the sense of the force between the black holes (see Ref. |P] and 
references therein). In the present paper, we restrict ourselves to the study 
of the system @. 

3 Auxiliary linear problem 

System of equations (|) is the compatibility condition for the following pair 
of matrix linear differential equations: 

Di-^ = U{uj,p,t)-^, U{uj,p,t) = ^ — , (10) 

p"^ — o;^ 

= V{uj, p, t)^, V{uj, p, t) = , (11) 

p"^ — cu^ 

Here Di and D2 are the commuting differential operators 

D, = dp-^^d^, D2 = dt-^^d^, (12) 
cu^ — p^ UJ^ — 

and is a complex parameter that does not depend on p,t. We also use the 
U — V-pair representation in which is a dependent parameter. To be more 
precise, let a; be a root of the equation 

cu^ + 2(t -k)LU + p^ = (13) 

where k is an independent parameter. Determined through (|T3|), uj satisfies 



a 2a;p 2cj2 

= ^ ^' = —2 (14) 

p^ — uj^ p^ — UJ^ 
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Passing from to \1/'(A;) = \E'(c<j(A;)), one obtains from (|T0 



n]) that 



dp^>' = u{uj{k),p,t)m\ dt^' = v{uj{k),p,t)m'. (15) 

We will denote the solution of (|T3|) by \E'(/c) missing the prime for brevity. It 
is worth mentioning that solving (|15D with the fixed branch of the root one 
finds the solution of ( p!0| , pj] ) only in the analyticity domain of uj{k). 

We will follow the general scheme for investigating integrable equations 
Ref. omitting many technical details. 

Eq. ( p!3D is invariant w. r. t. the transformation uj —>■ /oJ. Now suppose 
uj±{k) are root's branches such that \m.uj^ > and Imti;_ < 0. Then the cuts 
of u!± are half-lines (— oo,t — p] and [t + p, oo). Before proceed we mention 
some useful properties of uj±. Note that 

u4k) = -^ u+{k)=u_{k). (16) 

As p — > cxD the functions uj±{k,p,t) satisfy uniformly in k,t the following 
estimate 

uj±{k,p,t) = ±tp + {k-t) + 0{l/p). (17) 

Furthermore, introduce the monodromy matrixes T±{k, p, p'). By defini- 
tion, they are solutions to 

d,T^{p,p') = U{u;^{k),p,t)T^{p,p'), T^{p,p)=I. (18) 

Besides, zero caveture condition, 

Ut-V,+ [U,V]=0, 

reveals that 

dtTi{p,p') = V{u;^{k),p)T^{p,p')-T^{p,p')V{u;^{k),p'). (19) 

For symmetric real matrix g the system (0, |Tl]) remains invariant under 
the transformations 

^i/{u,p,t)-^g^-\^,p,t), ^{uj,p,t)^^{u,p,t) (20) 

UJ 
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where the tilde denotes transposition. Taking into account (|T6|) one gets 



T+{k,p,p')=g{p)Tz\k,p,p')g-\p'), T^k, p, p') = T^k, p, p') (21) 

Note that we regard the parameter as a parameter on the complex plane, 
not on the Riemann surface of the root, and understand equations ([T8| ) as 
two different ones. Let us remark also that T+(co',p, p') and T_{u, p, p') are 
solutions of (|T0|) as Imu > and Imu < respectively. 
Let us determine the Jost functions, 

v^±(A;,p) = jim T±(fc,p,p')e±(fc,p'), e±(A;,p') = J ^J^^^,^ ^ . (22) 

They are analytic in the k-plane with the cuts (— oo,t — p] and [t + p, oo). 
Moreover, from equation (|T9|) we derive that they are solutions of compatible 
system of equations (p^). Since limp_^oo e-(p)5'^^(p)e+(p) = / we see also that 

<il4k)=g^Z\k), ^4k) = ^4k). (23) 

It is worth mentioning that det '^±{k) = uj±{k). 

Recall that functions ^E'+(u;) and \E'_(ti;) are analytic solutions of ([TI1|), 
( pUD as Imuj > and as Imuj < respectively. Hence for real u there is a 
matrix T{k) such that 

2 2 

^+{u;) = ^.iu;)T{k), k = t+^^^. (24) 



Since ^±{uj) are coupled by transformation (poD and det\E'±(co') = lu, T{k) 
satisfies 

T{k) = f{k), f{k) = T-\k), detT{k) = 1, (25) 

by the other words it is a symmetric unitary matrix. It is interesting to note 
that if initial data {V — 1, X — 1, Vt, ^,t) are of compact support then T — I 
is also of compact support. In general case we have T{k) = / + 0(l/|/c|)as 
|A;| oo. 
Let 

X±(^) = ^±{u^)e-\u), e{u) = [ J M . (26) 
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It is possible to prove that x± are regular at = and X±(c>o) = /. The 
last condition and symmetry (20) reveal that 

^ = X±(0)( J p° ) (27) 

Since X-(0) = X+(0) from (|2|) we derive that Ti2{k) = T2i{k) = 0{l/\k\'^) 
at = oo. 

Boundary condition (|^ yields 

T/r M I f dtY\ 

UJ 
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Solving ([TT|) one obtains 



^±Mlp=o=(j T_(fc)=f+(fc) (28) 



Here T+ is analytic in k as ImA; > 0. Inserting (g^) in (^H) we come to 
identities: 

Tz\k)T+{k) = T{k) = f+{k)f-\k) (29) 
It follows from normalization of x± at c<j = oo that 

lim eiuj)TJk)e-\uj) = I (30) 

As a function of uj, T^{k{uj)) is analytic in domains: Imuj > 0, |ti;| > p 
and Imu; < 0, Ic^l < p. Then using (pij) and (120) one can easily check that 



^ i '^+i^)T+\k) lm^>0,|a;|>p . . 

\^4u)Tz\k) Ima; < 0, |cu| > p ^ ^ 

is analytic as |ci;| > p while 

(5 ,N / ^+(^)r-(A;) Im^ > 0, 1^1 < p , . 

- 1 v&-(^)r+(A:) Imc. < 0, < p ^^^^ 

is analytic as < p. On the circle = p, $± satisfy the conjugation 
condition, 

$+(cu) = $_(cj)G'(A;), ^(A;) = T+{k)f^{k) = T^{k)f+{k). (33) 
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Set x± = $±e ^. The matrices $+ and $_ are still coupled by transformation 
(pOl) . Using it and the fact that X-(oo) = J which follows from (^), we come 
to 

^? = X+(0)(J (34) 

Due to the boundary condition the limit, limp^oX+(0)) exists. On 
the other hand we have from (^) that 

limx+(0) = lim f J -^l^^]G{k)e-\u). 

The above limit exits if and only if Y{t) = Gi2(t)/G22{t)- Eventually, 

(nt) mt)]^(aii!> '■'■'■•"'"""1^1;^""""°"'" ^ ,35) 
I X{t) ) [ G.n(t) ) ^ ' 

where unimodularity of G{k) was used as well. 

We summarize the results of this section as follow. For any solution of the 
initial- value problem posed in Section 2 there exits analytic in upper fc-plane 
matrix, T+{k), that satisfies the conditions, 

T+T| = T|r+, ImT+T| = 

for real k. Here f is the Hermitian conjugate. It uniquely determines sym- 
metric unitary matrix, T{k), and symmetric real matrix, G{k). Then the 
inverse problem is reduced to the Riemann-Hilbert problem (RHP) ( p4D or 
(|5HD. In addition, there is one to one correspondence between G{k) and the 
data on the symmetry axis, for example Y{t) and V{t) (^). 



4 Einstein-Rosen waves 

In this section we apply the results of the previous section to a space-time 
with diagonal metric {W = 0). Set X = p^e"^'^. Then the system (H) reduces 
to one linear equation, viz 

- ^ + -TTT + = (36) 
dt^ dp^ pdp ^ ' 
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The above equation is the compatibihty condition of the following pair of 
scaler linear equations 



The solution of the system (^) is given by 

/ ge±(fc,p,t) Q \ 

with 

2p'^dpip 2p'uj±dtip 



e^{Kp,t) = -Jdp'l^^ + 



p'^ — — Uj. I 

The function k{uj) = t + maps the upper half-plane (lower half- 

plane) into fc-plane with the cuts (—00, t — p] and [t + p, 00). The function 
<jjj^{k, p'){uj_{k, p')) maps the fc-plane with the cuts (—00, t—p'] and [t+p', 00) 
into the upper half-plane (lower half-plane). Analyzing these mappings in the 
case when p < p' we conclude that 

uj+(k(uj + iO), p') - uj^{k{uj - iO), p') = 2i^J p'^ - {k ~ 
as |A; — t| < p' and zero otherwise. Note that |fc — t| > p for any real uo and 



uj±{k^ p') = k — t ± i^J p'"^ — {k — as p < \k — t\ < p' and real k. Therefore 
e+{k{uj + iO), p, t) - e_{k{uj - iO), p, t) = 2iA{k) (38) 
where _ 

,ik-t)dp^ + p'dtij 



A{k) = - J dp' 



\k-t\ 



The piece-wise analytic function 9{uj), 6{uj) = 0-{uj) as Imci; < and 
6{ijj) = 0+{uj) as Imu > 0, is a unique solution of the RHP (p^) with 9{oo) = 
0. This solution can be written as 

SM.ILaW (39) 

n J X — CO 
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Then we can represent the solution of (|3^) as 



V,(«.p) = i9(0) = ^/dA™ (40) 

— OO 

Changing the variable in the above integral one has 

m P) = - Ydk^M^ + 1 7 dk^^^= (41) 



5 Conclusions 

In this paper we outline the ISM scheme for the cylindrical symmetric waves 
model. However, to complete this scheme one needs to analyze the smooth- 
ness property of the matrices T{k) and G{k). We assume that for smooth 
initial data they are also smooth. The relation of the matrix G with the 
data on the symmetry axis allows to determine the behavior of the solution 
at the time infinities through the study of the asymptotic properties of G. 
The matrix T defines the behavior of the solution at the null infinities and 
we assume that it is uniquely restored from the data on the null infinities. 
An important open question is the derivation of the formula for the angle 
of deficit at the space infinity in terms of the matrix T or G. We plan to 
consider the problems mentioned in separate papers without promising to do 
it soon. 
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